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ABSTRACT 

We present a search for Trojan companions to 25 transiting exoplanets. We use the technique of Ford & 
Gaudi, in which a difference is sought between the observed transit time and the transit time that is calculated 
by fitting a two-body Keplerian orbit to the radial-velocity data. This technique is sensitive to the imbalance 
of mass at the L4/L5 points of the planet-star orbit. No companions were detected above 2a confidence. The 
median 2cr upper limit is 56 M®, and the most constraining limit is 2.8 M ffi for the case of GJ 436. A similar 
survey using forthcoming data from the Kepler satellite mission, along with the radial-velocity data that will 
be needed to confirm transit candidates, will be sensitive to 10-50 Mq Trojan companions in the habitable 
zones of their parent stars. As a by-product of this study, we present empirical constraints on the eccentricities 
of the planetary orbits, including those which have previously been assumed to be circular. The limits on 
eccentricity are of interest for investigations of tidal circularization and for bounding possible systematic errors 
in the measured planetary radii and the predicted times of secondary eclipses. 

Subject headings: techniques: transit photometry, radial velocities — extra-solar planets, trojans — eccentricity 



1. INTRODUCTION 

Trojan companions are bodies in a 1:1 mean-motion reso- 
nance with a planet, librating around one of the two triangular 
Lagrange points (L4 and L5) of the planet's orbit around the 
star. The archetypal example is the population of Trojan as- 
teroids in resonance with Jupiter. Trojan companions to Nep- 
tune and Mars have also been detected (Sheppard and Tru- 
jillo 2006, Rivkin et al. 2007). Another interesting example is 
the pair of Saturnian satellites Calypso and Telesto, which are 
in 1 : 1 resonance with their fellow satellite Tethys (Reitsema 
1981). The presence of Trojan companions and their orbital 
and physical characteristics have been considered as clues to 
processes in planet formation and migration. Several recent 
studies have examined the capture and survival of Trojans in 
the context of suspected changes in the orbital architecture of 
the Solar system (Morbidelli et al. 2005, Chiang and Lithwick 
2005, Kortenkamp et al. 2004). 

Although the Trojan-to-planet mass ratios in the Solar sys- 
tem are very small (mj jmp ~ 10~ 7 for Jupiter), it is conceiv- 
able that Trojans with much higher mass ratios exist in exo- 
planetary systems. For circular orbits, even very massive Tro- 
jans can be dynamically stable. Laughlin & Chambers (2002) 
explored the viability of Trojans with mass ratios of unity (i.e., 
co-orbital planets of equal mass), finding that such configura- 
tions can be dynamically stable over time scales comparable 
to or longer than stellar lifetimes. More generally, the sta- 
bility of the L4/L5 points depends on the orbital eccentricity 
and the relative masses of the Trojan, planet, and star (see, 
e.g., Nauenberg 2002, Dvorak et al. 2004). For many of the 
known exoplanets, considerations of dynamical stability al- 
low for massive Trojan companions. For example, at least 
7 of the known gas giant planets that are within the habit- 
able zones of their parent stars could have dynamically-stable, 
terrestrial-mass Trojan companions (Schwarz et al. 2007). 

Several methods have been proposed to detect Trojan com- 
panions to exoplanets. A Trojan may be massive enough to 
perturb the stellar motion by an amount that is detectable in 
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the radial-velocity (RV) orbit of the star (Laughlin & Cham- 
bers 2002). A Trojan in a nearly edge-on orbit may be large 
enough for its transit to be detected photometrically (see, e.g., 
Croll et al. 2007). For a transiting planet, the gravitational per- 
turbations from a Trojan companion may cause a detectable 
pattern in the recorded transit times (Ford & Holman 2007). 
Alternatively, Ford & Gaudi (2006) proposed comparing the 
measured transit times with the times that would be expected 
based only on the RV data and the assumption of a two-body 
orbit. 

An important virtue of the latter technique is that a sensitive 
search for Trojans can be performed using only the RV and 
photometric data that are routinely obtained while confirm- 
ing transit candidates and characterizing the planets. This is 
in contrast to the first three methods, for which it is generally 
necessary to gather new and highly specialized data (very pre- 
cise RVs, continuous space-borne photometry, and a long se- 
quence of precisely-measured transit times, respectively). For 
example, Ford & Gaudi (2006) and Narita et al. (2007) placed 
upper limits on Trojan companions of approximately Neptune 
mass to the transiting planets HD 209458b, HD 149026b and 
TrES-lb, using data gathered for other purposes. 

In this paper, we present a search for Trojan companions 
to 25 known transiting exoplanetary systems for which suit- 
able data are available, using the method of Ford & Gaudi 
(2006, hereafter, "FG"). This paper is organized as follows. 
The method is described in § 2. The compilation and analysis 
of the data is described in § 3. The results are given in § 4. 
These results are summarized and discussed in § 5, which also 
looks ahead to the prospects for a similar search using data 
from the Kepler mission (Borucki et al. 2008). 

As will be explained in § 2, the orbital eccentricity of the 
planet-star orbit affects the interpretation of the data. Hence, 
a necessary part of our analysis was the determination of the 
orbital eccentricity for each system, or the justification of the 
common assumption that the orbit is circular due to tidal ef- 
fects. These issues are investigated systematically in § 3. Our 
findings may be of interest independently of our results on 
Trojan companions, not only because of the connection to the 
theory of tidal circularization, but also because the orbital ec- 
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centricity affects estimates of the planetary radius via transit 
photometry, as well as the predicted times of planetary occul- 
tations (secondary eclipses). We discuss these points in § 5. 

2. METHOD 

The basic idea of the FG method is to compare the mea- 
sured transit time with the expected transit time that is cal- 
culated by fitting a two-body Keplerian orbit to the RV data. 
We will denote by to the observed transit time, and by tc the 
calculated transit time, in which the calculation is based on 
fitting a two-body Keplerian orbit to the RV data. The pres- 
ence of a Trojan companion as a third body would cause a 
timing offset At = to~tc- 

This is most easily understood for the case of a planet on a 
circular orbit. In such a case, if there is no Trojan companion, 
the force vector on the star points directly at the planet, and 
the observed transit time to coincides with the time ty when 
the orbital velocity of the star is in the plane of the sky (i.e., the 
time corresponding to the null in the RV variation). If instead 
there is a single Trojan located at L4 or L5 (or librating with 
a small amplitude), then the force vector on the star does not 
point directly at the planet; it is displaced in angle toward 
the Trojan companion, given by tan(0) ~ s/3e/(2 — e) where, 
e = mp /(mp+mp) for a Trojan mass mp and a planet mass nip 
(Ford & Gaudi 2006). As a result, to occurs earlier or later 
thanfy, and the time difference is given by At = ±4>P/2ir. For 
small values of the Trojan-to-planet mass ratio, the magnitude 
of to— ty is proportional to mp, (Ford & Gaudi 2006): 
P \ ( m T \ f0.5M Jap s 



At ~ ±37.5 min 



(1) 



mp = mp 



(2) 



3 days J \ lOM^ J \mp + mp 
The positive sign corresponds to a mass excess at the L4 point 
(leading the planet) while the negative sign corresponds to a 
mass excess at the L5 point (lagging the planet). Thus, given a 
At, the mass excess can be estimated using Eq. (Q3, assuming 
small Trojan-to-planet mass ratio. 
More generally, the mass excess is given by: 
2tan(27rAf/F) \ 
,V3-|tan(27rAf/F)|/ 
For an eccentric two-body orbit, the transit time does not 
generally coincide with the time of null RV variation, and 
hence in general tc^ty. To first order, tc~ty « (e cosw)P/27r, 
where e is the eccentricity and ui is the argument of peri- 
center, and hence one may use the statistic At = to~ty - 
(e cos uj)P /2ir to search for Trojan companions. This is how 
the problem was described by FG, although we find it use- 
ful to cast the problem more generally as a comparison be- 
tween to and tc- We emphasize here that to depends solely on 
photometric observations of transits, while tc depends almost 
entirely on RV observations. 2 

Specifically, one calculates tc by fitting a two-body Keple- 
rian orbit to the RV data and calculating the expected transit 
time based on the the fitted orbital parameters (see, e.g., Kane 
et al. 2008). The true anomaly (/) corresponding to the transit 
time is 

7T 

(3) 
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from which the eccentric anomaly E can be calculated using 



E 

tan — = 

2 



l-e f 
1+e 2 
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2 As explained in § [3] the only sense in which tc depends on photometric 
data is that we used the photometrically-determined orbital period P when 
fitting the RV data, to reduce the number of free parameters. 



which in turn leads to the mean anomaly M of the transit using 
Kepler's equation, 



M = E-esinE. 



(5) 



Finally, the calculated transit time tc is obtained from the def- 
inition of the mean anomaly, M = 2ir(t-tp)/P, where tp is the 
time of pericenter passage. 

Our basic procedure is therefore to determine t Q from pub- 
lished transit ephemerides, calculate t c by fitting a two-body 
Keplerian orbit to the available RV data, and calculate Af = 
to~tc- For circular orbits, we use Eq. (f2]i to determine the 
Trojan companion mass mp corresponding to a given value of 
At. For eccentric orbits, the relationship between tc and mj 
is determined using direct numerical integrations of 3-body 
systems using a Bulirsch-Stoer algorithm (Varadi et al. 1996). 
In these integratons, we hold fixed P, e, u, and the stellar 
mass m$ at the values given in the literature, and select a Tro- 
jan mass mj and planetary mass mp such that the RV semi- 
amplitude (Nauenberg 2002) 



K = 



2ttG 



1/3 



y mp + mj + mpmp 
(ms+mp + m.T) 2 / 3 V l—e 2 
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is equal to the observed value. Hence we simulate the case 
in which the observed RV variation is due to the combined 
force of a planet and a Trojan, rather than a planet alone, but 
the RV data alone are insufficiently precise to discern the dif- 
ference. 3 We compute the transit time t c , repeat the analysis 
for an increasing sequence of mp, and fit a polynomial func- 
tion to the resulting relationship tc(mp)-tc(Q). We found a 
quadratic function, mp = a\ At + «2(Af) 2 , to give a good fit to 
the results. Taking mp to be in Earth masses and Af in min- 
utes, the coefficients (a x ,ai) are (0.044,-1. 17x 10" 5 ) for GJ 
436b and (6.787,-0.001) for XO-3b. For the cases of HAT-P- 
2b and HD 17156b, we find that even very low-mass Trojan 
companions are dynamically unstable, owing to the large or- 
bital eccentricities (see § 14.2b . Thus, for those systems, the 
requirement of dynamical stability is more constraining than 
the empirical upper limit on mp based on the FG method. (As 
will be described in § 4, this also proved to be true for XO-3b 
based on the current data.) 

3. DATA ANALYSIS 

The RV data were taken from the available literature on 
each system. The references are given in Table 1. These 
data were generally obtained for the purpose of discovering or 
confirming the planet, although in a few cases the data were 
obtained for other reasons, such as precisely measuring the 
orbital eccentricity (Laughlin et al. 2005) or for measuring 
the Rossiter-McLaughlin effect (Winn et al. 2006). Regard- 
ing the latter, the data that were obtained while a transit was 
in progress were not used, to avoid the needless complication 
of incorporating the Rossiter-McLaughlin effect into the RV 
model. However, in those cases the investigators usually gath- 
ered additional data outside of the transit which are useful for 
refining the spectroscopic orbit. 

Our RV model for an eccentric Keplerian orbit has \+N free 
parameters, where is the number of independent data sets. 
Those parameters are the projected planet mass (mpsinO, or- 
bital eccentricity (e), argument of pericenter {uS), calculated 
time of midtransit (tc), and a constant additive velocity (7) 

3 We verified that this discernment is indeed impossible, for mp /nip < 0.5, 
for the systems with eccentric orbits considered in this paper. 
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TABLE 1 
Description of Data 



System 


N, [ ' 


Jitter [m s ] 


<x v [m s" 1 ] 


K[m s~'] 


C — 1 / £77 ™ m / ™ i 
s / \tn j j mp ) 


References 


HD 209458 


55 


1.4 


4.9 


83.3 


77.7 


1,2 


HD 17156 


24,8 


3.3,3.4 


3.7,6.2 


273.4 


72.0 


47,48,49,50 


HAT-P-7 


8 


6.5 


6.7 


213.6 


55.3 


51 


HD 189733 


16,44 


0.0,0.0 


12.0,3.0 


201.3 


40.6 


4 


TrES-2 


11,5 


0.0,0.0 


7.5,5.8 


181.5 


32.2 


18,19 


HAT-P-3 


9 


5.1 


5.8 


98.7 


31.3 


10 


HAT-P-4 


9 


4.4 


5.0 


80.8 


29.6 


1 1 


WASP-3 


6 


0.0 


13.8 


247.7 


26.9 


27 


WASP-5 


1 1 


11.6 


21.8 


276.4 


25.8 


29 


TrES-3 


1 1 


30.0 


31.6 


370.4 


23.8 


20 


HAT-P-6 


13 


10.1 


11.2 


116.2 


23.0 


13 


GJ436 


52 


3.4 


4.2 


18.3 


19.4 


45,46 


WASP-4 


13 


22.3 


28.0 


240.3 


18.9 


28 


HD 149026 


16 


5.4 


6.1 


46.4 


18.7 


3 


HAT-P-2 


13,10,7 


34.8,88.7,21.4 


35.5,104.5,24.7 


980.0 


16.8 


7,8,9 


WASP-2 


7 


16.6 


17.7 


157.6 


14.4 


25,26 


WASP-1 


7,5 


2.7,0.0 


3.9,13.3 


127.8 


12.8 


22,23,24 


XO-3 


10,10 


0.0,0.0 


171.3,159.0 


1486.2 


12.3 


34 


HAT-P-1 


15,8 


4.4,0.0 


6.4,7.4 


59.0 


11.9 


5.6 


TrES-4 


4 


0.0 


10.8 


98.3 


11.1 


21 


TrES-1 


7,8,5 


0.0,0.0,0.0 


12.1,14.6,3.4 


112.0 


9.7 


14,15,16,17 


CoRoT-Exo-1 


9 


34.0 


47.5 


190.9 


7.4 


52 


CoRoT-Exo-2 


8,4,3,9 


60.0,60.0,0.0,0.0 


68.4,61.7,27.7,19.0 


594.4 


7.3 


53,54 


XO-2 


9 


15.6 


25.0 


84.1 


6.2 


33 


HAT-P-5 


8 


33.4 


37.7 


134.0 


6.2 


12 


OGLE-TR-182 


20 


29.7 


59.5 


120.0 


5.5 


43 


OGLE-TR-113 


8 


83.5 


93.1 


286.1 


5.3 


41 


OGLE-TR-211 


20 


23.7 


55.4 


82.0 


4.1 


44 


OGLE-TR-56 


11 


89.0 


153.1 


268.3 


3.6 


37,38 


OGLE-TR-111 


8 


0.0 


40.2 


78.0 


3.4 


39,40 


OGLE-TR-132 


5 


51.0 


68.7 


167.0 


3.3 


42 


OGLE-TR-10 


9 


0.0 


63.2 


80.0 


2.3 


35,36 


XO-1 


4.6 


0.0,0.0 


65.1,16.8 


120.1 


2.2 


30,31,32 



REFERENCES. — (1) Laughlin et al. 2005a; (2) Winn et al. 2005; (3) Wolf et al. 2007; (4) Winn et al. 2007a; (5) Bakos et al. 2007a; (6) Winn et al. 2007b; (7) Bakos 
et al. 2007b; (8) Winn et al. 2007c; (9) Loeillet et al. 2007; (10) Torres et al. 2007a; (11) Kovacs et al. 2007; (12) Bakos et al. 2007c; (13) Noyes et al. 2008; (14) Alonso et 
al. 2004; (15) Laughlin et al. 2005b; (16) Narita et al. 2007a; (17) Winn et al. 2007d;(18) O'Donovan et al. 2007a; (19) Holman et al. 2007; (20) O'Donovan et al. 2007b; (21) Man- 
dushev et al. 2007; (22) Collier Cameron et al. 2007; (23) Stempels et al. 2007; (24) Charbonneau et al. 2007; (25) Winn et al. 2008; (26) Charbonneau et al. 2007; (27) Pollaco et 
al. 2007; (28) Wilson et al. 2008; (29) Anderson et al. 2008; (30) McCullough et al. 2006; (31) Holman et al. 2006; (32) Wilson et al. 2006; (33) Burke et al. 2007; (34) Johns-Krull et 
al. 2007; (35) Konacki et al. 2005; (36) Pont et al. 2007a; (37) Torres et al. 2004; (38) Pont et al. 2007b; (39) Pont et al. 2007c; (40) Winn et al. 2007e; (41) Bouchy et al. 2004; (42) Bouchy 
et al. 2004; (43) Pont et al. 2007d; (44) Udalski et al. 2008; (45) Maness et al. 2007; (46) Gillon et al. 2007; (47) Fisher et al. 2007; (48) Narita et al. 2007b; (49) Gillon et al. 2007; (50) Ir- 
win et al. 2008; (51) Pal et al. 2008; (52) Barge et al. 2008; (53) Alonso et al. 2008; (54) Bouchy et al. 2008. 
d Multiple values represent multiple data sets available for the system. 



for each data set. In practice we use parameters ecosw and 
esinw instead of e and uj because for small e, the errors in 
ecosw andesinoj are uncorrected (see, e.g., Winnetal. 2005, 
Shen & Turner 2008). The orbital period P is held fixed at 
the photometrically determined value, but of course the tran- 
sit time tc is not constrained by the photometric data, since 
it is the difference between tc and the actual transit time to 
that we are trying to measure. In two cases for which a long- 
term acceleration has been identified in the RV data (GJ 436b 
and CoRoT-Exo-lb), we include an additional free parameter, 
7 4 We assign a different 7 to each RV data set, to allow for 
telescope-specific velocity offsets. The stellar masses for the 
systems were taken from the homogeneous analysis of Torres 
et al. (2008) when possible, and otherwise from the discovery 
paper. 

3.1. Estimation of Jitter 

4 For the remaining systems, we have assumed that the acceleration term 
is zero. Any real acceleration (and any other failures of the single-Keplerian 
model) will appear as "noise" in our a nalysis and will be reflected in a larger 
estimate of stellar "jitter" (see § |3.1> . We find that our results for At do 
not depend much on whether or not a long-term acceleration is allowed as 
an additional parameter, because the error in this parameter is not strongly 
correlated with the error in Af. 



For each system, we first analyzed the data with the goal 
of reproducing the quoted results in the literature. We fitted a 
Keplerian model to the RV data by minimizing the \ 2 statis- 
tic using the AMOEBA algorithm (Press et al. 1992). The 
initial conditions for the free parameters were taken to be the 
literature values, and for consistency, for this step we used the 
exact same choices of P and nis as in the literature. We define 
X 2 as 

x 2 = E(^) 2 , CD 

where v,- o and v, c are the observed and calculated radial ve- 
locities, respectively, and 07 is the corresponding uncertainty. 
The uncertainty should include the statistical uncertainty cr sta t, 
as well as the systematic error cT sys due to unmodeled instru- 
mental systematic errors and intrinsic variations of the stellar 
photosphere, often referred to as "stellar jitter" in this con- 
text (Wright 2006). To estimate the appropriate values of 
o~i for this project, we determined the value of cr sys such that 
X 2 /Ndof= 1 when using 



<Ti= V^stat+^ys (8) 
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in Eq. (jTj. Our estimates of the stellar jitter using this proce- 
dure are given in Table 1 . 

3.2. Data Selection 

Our desire was to perform as wide a search as possible, us- 
ing all publicly available data, but for some systems the RV 
data is so sparse that meaningful constraints cannot yet be ob- 
tained. To guide our selection of systems, we used a figure-of- 
merit based on the results of the Fisher information analysis 
presented by FG. Those authors showed that in the limit of 
continuous RV sampling with uniform errors, the uncertainty 
in Af will approach a A , = {l/2-n: 2 N,) l/2 P(j v /K, where N v is 
the number of radial-velocity data points and er v is the error 
per point. For a circular orbit the corresponding uncertainty 
in mj will approach a mj = {%/3N v ) l l 2 mpa v /K. Thus, it is 
possible to estimate the expected uncertainty in At and nip, 
given the published system parameters, without fitting the ac- 
tual data. We calculated the figure of merit 



S = Z ~ L / a (m T /m P ) 



CJy 



(9) 



for each system in the literature at the outset of this project, 
and ranked the systems accordingly. Table 1 shows £ for all 
the systems, along with N v , o~ v , and K. For systems where 

1 /2 

multiple data sets are available, the effective a v /N v is ob- 
tained by adding in quadrature the corresponding terms from 
the different data sets. A higher value of £ reflects better qual- 
ity of data. From MCMC analyses of all the systems, we find 
that for systems with £ < 6, the fitting algorithm is suscepti- 
ble to poor convergence and allows for unphysical parameter 
ranges. Since the scientific return on such low-£ systems is 
comparatively poor we decided to remove them from consid- 
eration rather than eliminate these fitting problems. In what 
follows we focus exclusively on the 25 systems with £ > 6. 

3.3. Assumptions for orbital eccentricity 

As explained in § 2, the orbital eccentricity affects the cal- 
culation of tc and also affects the relation between Af and the 
possible Trojan companion mass mj. Hence it is imperative 
to consider the possiblity of eccentric orbits. Most of the cur- 
rently known transiting planets are in very close orbits, where 
the effects of tidal interactions between the star and planet — 
and orbital circularization in particular — are expected to be 
significant (Rasio et al. 1996, Trilling 2000, Dobbs-Dixon et 
al. 2004). A common practice is to assume that, in the ab- 
sence of positive evidence for an eccentric orbit, the orbital 
eccentricity has been reduced to insignificance by the action 
of tides. 

If the assumption of a circular orbit could be justified, it 
would be advantageous for the present study because it would 
remove 2 free parameters from the Keplerian model (e and to) 
and thereby strengthen the determination of the other param- 
eters, including the key parameter tc- Our approach was to 
assume the orbit to be circular only when (1) a circular orbit 
is consistent with the RV data, (2) the estimated stellar age 
is more than 20 times larger than the estimated timescale for 
tidal circularization, and (3) no constraint on ecosui is avail- 
able because no planetary occultations (secondary eclipses) 
have been observed. These points are explained in detail in 
the paragraphs to follow. 

To test whether the RV data are consistent with a circu- 
lar orbit, we fitted a Keplerian model to the RV data using 



a Markov Chain Monte Carlo (MCMC) technique, employ- 
ing a Metropolis-Hastings algorithm within the Gibbs sam- 
pler (see, e.g., Tegmark et al. 2004; Ford 2005; Holman et 
al. 2006; Winn et al. 2007a). For this step, the free parameters 
were nip sin i, ecosw, esincj and a 7 for each data set. Uni- 
form priors were used for all parameters. The fitting statistic, 
X 2 , was defined in Eq. ©. A single chain of ~ 10 6 links 
was used for each system. The jump sizes for the various 
parameters were set such that the acceptance rate for each pa- 
rameter was ~20%. For each parameter, we found the mode 
of the a posteriori distribution (marginalized over all other 
parameters), and the 68.3% confidence interval, defined as 
the range that excludes 15.9% of the probability at each ex- 
treme of the a posteriori distribution. For cases when ecosui 
and esinu! were both consistent with zero, we also found the 
95.4%-confidence upper limit on e. Table 2 gives the results. 
All of the systems were found to be consistent with a circular 
orbit, except for the 4 well-known eccentric systems GJ 436, 
HAT-P-2, HD 17156, and XO-3. 

For the estimated timescale for tidal circularization, we 
used (Goldreich & Soter 1966): 



4 J a 3 
Tcirc " 63 G Um s 



1/2 



nip / a 
m s \R„ 



(10) 



which is based on the highly simplified, widely-used model 
of tidal dissipation in which the tidal bulge experiences a con- 
stant phase lag due to tidal friction. Here, a is the orbital sep- 
aration, and R p is the radius of the planet. The dimensionless 
number Q is inversely proportional to the dissipation rate. In 
the solar system, Jupiter is thought to have Q ~ 10 5 (Ioannou 
& Lindzen 1993) to the extent that this simplified model is ap- 
plicable. For our purpose, a necessary condition for assuming 
the orbit to be circular was that r^/r c i rc > 20, i.e., there have 
been at least 20 e-foldings of tidal circularization, according 
to this model. In calculating T c ; rc we assumed Q = 10 6 , which 
is conservative in the sense that a larger Q corresponds to a 
longer calculated timescale for circularization, and a smaller 
risk that we are assuming a circular orbit when this assump- 
tion is not justified. 

For a few systems, additional constraints on the orbital ec- 
centricity are available because a planetary occultation (sec- 
ondary eclipse) has been observed. For small eccentricities, 
the variables ecosw and esinw are directly related to the in- 
terval between the transit and occultation, and the relative du- 
rations of those two events (Kallrath & Milone 1999, Char- 
bonneau et al. 2005): 



= ^ 1 focc" 



ftra — ~ 



6^-00 

e, ra +e 



(ii) 



(12) 



where f tra and ? occ are the times of transit and occultation, and 
0tr a and OCC are the corresponding durations. In all cases 
to date, the bounds on esinw that follow from this relation 
are weaker than bounds from the RV data (see, e.g., Winn 
et al. 2005), and the bounds on ecosw are more constrain- 
ing. Thus, for those cases in which an occultation has been 
observed, we add a term to our \ 2 statistic to enforce the cor- 
responding constraint on ecosw: 



N 2 

^ 2 „y^ f vo-vc \ | / (ecosuj)o-(ecosuj)c s 

1 \ / \ ®e cosuj , 

n=l N / x / 



(13) 







Trojan Companions to Transiting Exoplanets 




5 






TABLE 2 

Inferred orbital eccentricities and related parameters 






System 


ecostj 




' cue 


T*/r circ 



CoRoT-Exo-1 


+0. 


CoRoT-Exo-2 


-0. 


GJ436 


+0. 


HAT-P-1 


+0. 


HAT-P-2 


-0. 


HAT-P-3 


+0 


HAT-P-4 


-0. 


HAT-P-5 


+0. 


HAT-P-6 


+0. 


HAT-P-7 


-0. 


HD 149026 


-0. 


HD 17156 


-0. 


HD 189733 


+0. 


HD 209458 


+0. 


OGLE-TR-10 


+0. 


OGLE-TR-111 


+0. 


OGLE-TR-113 


-0. 


OGLE-TR-132 


+0. 


OGLE-TR-182 


-0. 


OGLE-TR-211 


+0. 


OGLE-TR-56 


+0. 


TRES-1 


+0. 


TRES-2 


+0. 


TRES-3 


+0. 


TRES-4 


+0. 


WASP-1 


+0 


WASP-2 


-0. 


WASP-3 


-0. 


WASP-4 


+0. 


WASP-5 


+0. 


XO-1 


+0. 


XO-2 


+0. 


XO-3 


+0. 



,+0.038 

-0.025 
1+0.006 

-0.006 
q+0.021 

-0.029 
(.+0.005 
1 -0.006 
3+0.053 



0.013 



0.044 



-0.014 
c+0.095 
l '-0.()95 
3+0.016 

Urn 



J -0.886 
3+0.616 
-0.028 
+0 092 

-Ml 
'M 
+8:fi 

-0.170 
3+0.279 
3 -0.790 
3+0 002 
J -0.002 
-1+0.015 

-fi.QII 



-0 642 
fr+0.031 

-0.038 
1+0.331 

M 



--0 017 
7+0.064 



-0.073+0;" 
+0 054 + " : "-5 

™- UJ -0.027 

40.00**;°* 
-0 059+" 1114 

+0 033+ 062 
u -0.103 

-0 n54 +0 '° 5 * 

0.039 q 

+0 042+° : ° 34 
-0.034 

+0 000 +iiBI<5 

+n i o+0.042 
-t-u.iuy 06g 

+n <: 7 3+0.006 
tu. J / ->_ .006 

+0.008+; ' 
+0 436 +ii ' 3 " 

+U ' 4JD -(I 264 

+0 099+° 072 

™' u ' -0 559 

+0 1S2+ ' 104 
+0.279+ - 084 
40.352^1 

4-0 144+0244 
+U.lf-t_ H 244 
+ A 51O+0.25 1 

-0 039+ 030 
-0 024 +!i < il§ 
-0 031 +B M6 

u - u;u -0.040 
+0 4^4-0-054 

+U.tJt_ H 590 

+0 009+° 035 
+u.uuy 039 

-0.016+ofs 

+0.015^| 

+ 0.026 + °;°| 
40.066 + | 
-0 216 +ii:i6 ° 
-0.063+™ 4 



< 0.284 

< 0.101 

°- 138 ±o:oo7 

< 0.099 
520+ 004 

< 0.194 

< 0.123 

< 0.442 

< 0.101 

< 0.038 

< 0.179 
fi69 +0.008 

< 0.024 

< 0.028 

< 1.000 

< 0.964 

< 0.417 

< 0.993 

< 0.960 

< 0.858 

< 0.998 

< 0.084 

< 0.078 

< 0.101 

< 0.859 

< 0.088 

< 0.547 

< 0.098 

< 0.096 

< 0.088 

< 0.290 

< 0.515 

n 77Q+O.OI6 
u - zzy -0.018 



8.00 
0.50 
6.00 
2.70 
2.60 
1.50 
4.60 
2.60 
2.30 
2.20 
1.90 
5.70 
6.80 
3.10 
3.20 
8.80 
13.20 
1.20 
2.00 
2.00 
3.20 
3.70 
5.00 
0.60 
2.90 
3.00 
5.60 
2.10 
2.00 
2.00 
1.00 
5.80 
2.70 



0.00 
0.02 
1.15 
0.37 
57.59 
0.31 
0.09 
0.09 
0.32 
0.05 
1.10 
2152.22 
0.05 
0.11 
0.13 
0.51 
0.01 
0.02 
0.83 
0.15 
0.00 
0.19 
0.07 
0.00 
0.05 
0.02 
0.05 
0.02 
0.00 
0.03 
0.43 
0.12 
2.98 



3179.24 
29.86 
5.23 
7.30 
0.05 
4.78 
50.18 
27.98 
7.26 
47.88 
1.73 
0.00 
136.62 
27.22 
25.54 
17.39 

1204.61 
72.76 
2.41 
12.91 

1139.56 
19.75 
75.20 
137.03 
54.60 
129.19 
116.01 
93.68 

931.94 
77.79 
2.32 
48.43 
0.91 



a 95.4 % confidence limits on eccentricity. For four systems which are clearly eccentric, the mode and the 68.3% confidence limits are reported. 

k Nominal age of the system in Gyr : Taken from Torres et al. 2008, when possible, and from the discovery papers for systems not analyzed by Torres et al. 2008. For Corot-Exo-1, we 
assume a nominal age of 8 Gyr based on the reasoning that is has to be a fairly old main-sequence star (see Barge et al. 2008). 
c Circularization time-scale (in Gyr) for the planetary orbit, assuming Q = 10 6 (see the text). 



TABLE 3 

Constraints on ecosu from Secondary Eclipse Observations 



System 


ecosLo" 


Reference b 




HD 209458 


0.0002 ±0.0021 


Deming et al. 


2005 


HD 189733 


0.0012 ±0.0002 


Knutson et al. 


2007 


HD 149026 


-0.0011 ±0.0009 


Harrington et al. 


2007 


TrES-1 


0.0030 ±0.0019 


Charbonneau et al. 


2005 


GJ436 


0.1 346 ±0.0059 


Deming et al. 


2007 



a Evaluated from the transit ephemeris and the secondary eclipse 
time using Eqn. il lb . 

b Reference in literature from which the secondary eclipse time was 
obtained. 



where (ecosw) and (ecosw) c are the "observed" and calcu- 
lated values of ecosuj. By "observed" we mean the value that 
follows from the measured interval between transits and oc- 
cultations when inserted into Eq. ( fTTT i. The 5 systems in this 
category, and the constraints on ecosw, are listed in Table 3. 
In those cases, because such a powerful empirical constraint is 
available, we do not assume the orbit to be circular even if the 
RV data are consistent with a circular orbit and t*/ T c ; rc > 20. 

3.4. Constraints on the masses of Trojan companions 



For each system, after determining the appropriate level of 
stellar jitter and deciding whether or not the assumption of a 
circular orbit is justified, we determined the key parameter tc 
and its uncertainty using the same MCMC code that was de- 
scribed in the previous section. In all cases the free parameters 
included nip sin i, 7, and tc- In cases for which the a circular 
orbit was not assumed, we also fitted for ecosw and esinw. 
For the special cases of GJ 436 and CoRoT-Exo-1, a velocity 
gradient 7 was also included as a free parameter. The basic 
fitting statistic, % 2 , was defined in Eq. Q, and for the systems 
in Table 3 an a priori constraint on e cos ui was applied as in 
Eq. dTill- 
To determine the photometric transit time t , we used the 
most precise published photometric ephemeris to compute a 
predicted transit time close to the midpoint of the RV time 
series. We then computed the key parameter At = to- tc, 
the difference between the photometrically observed transit 
time and the transit time calculated from the RV data assum- 
ing zero Trojan mass. In all cases, the uncertainty in to is 
negligible in comparison to the uncertainty in t c - The results 
for At are translated into constraints on the Trojan mass m T 
using Eq. (|2]i for circular orbits, and using the numerical inte- 
grations described in §|2]for eccentric orbits. 
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the result for Af was consistent with zero within 2a. The 
system that was closest to a 2a detection was WASP-2, for 
which Af = -123±g minutes. The result for WASP-2 is there- 
fore worth following up with additional RV data. However, 
in a sample of 25 systems, even if Af is always consistent 
with zero, one expects approximately one 2a outlier. Hence 
our survey has not produced compelling evidence for a Trojan 
companion in this ensemble. 

The 2a upper limits on nij and on nij/mp are shown in 
Figures [3] and |4] respectively. The systems are ordered from 
least-constrained to best-constrained, going from left to right. 
The median upper limit on mj is 56 M ffi , with the most con- 
straining limit of 2.8 M ffi holding for the Neptune-sized planet 
GJ 436. Such a powerful upper limit is possible in this case 
because of the small stellar and planetary masses, and the co- 
pious RV data that is available for this system. The median 
upper limit on the mass ratio mj /nip is 0. 1 . 

It is possible to compare our results to those obtained pre- 
viously for 3 particular systems. For HD 209458, FG found 
Af = 13 ± 9 minutes and we find 2 + _\ y min. For HD 149026, 
FG found Af = -13 ± 27 min and we find 26^g min. For 



FIG. 1. — Constraints on ecosu and esmui based on fitting a Keplerian 
orbit to the RV data, and using constraints from the observed time of the 
secondary eclipse when available. The four clearly eccentric systems are 
labeled. The dotted circles show the eccentricity contours corresponding to 
e = 0.2, 0.4 and 0.6. 



-0.6 -0.4 -0.2 -0.0 0.2 0.4 0.6 TrES-1, Narita et al. (2007) found Af = -3.2 ± 11.8 min, as- 
e cosfw) suming a circular orbit, and we find -4^[ j min, allowing the 

v ' orbit to be eccentric but using the constraint on e cos u from 

secondary eclipse. These results are all consistent with zero 
with approximately the same range of uncertainty. Minor dif- 
ferences in the quoted central values are probably attributable 
to minor differences in the fitting procedures and in report- 
ing median values of the a posteriori distributions rather than 
modes. We also find our uncertainties to be in general agree- 
ment with the forecasted uncertainties based on the Fisher in- 
formation analysis of FG. 

4.2. Considerations of dynamical stability 

For a planet on a circular orbit, non-librating Trojan com- 
panions are stable as long as the masses satisfy the condition 
(Laughlin & Chambers 2002): 




mp+niT 
{nis + nip + niT) 



< 0.03812, 



(14) 



where m$, nip and nij are the masses of the star, planet and 
Trojan companion, respectively. This criterion allows for Tro- 
jan "companions" that are just as massive as the planet itself, 
even for planets as massive as 10 Mj up around a Sun-like star. 

However, the condition for stability of Trojan companions 
depends strongly on the eccentricity of the orbit. Nauenberg 
(2002) reported just such a study, showing the stability do- 
main of bodies in 1 : 1 resonance as a function the eccentricity 
of the orbit and the Routh parameter (7«) given by: 



1R ] 



ntstrip + nipnip + mstrvp 



(15) 



-0.10 -0.05 0.00 0.05 
e cos(w) 



0.10 



FIG. 2. — The inner region of FigfTJ The dotted circles show the eccentricity 
contours corresponding to e = 0.05 and 0.1. 



4. RESULTS 

4.1. Constraints on Trojan Masses 

Table 4 gives the 68.3% (ler) confidence intervals for At 
and nij for all 25 systems under consideration, as well as the 
95.4% (2a) upper limits on nij and nij /nip. In all the cases, 



(ms+mp + niT) 2 

Given the masses of the three bodies, and the eccentricity of 
the system, one can calculate jp and determine from Fig. 5 
of Nauenberg (2002) whether or not the system is stable in 
1 : 1 resonance. Conversely, given the eccentricity of the sys- 
tem, Fig. 5 of Nauenberg (2002) gives the the maximum 7^ 
allowed for stability which, along with nig and nip, gives the 
maximum Trojan mass allowed in the system. 

We calculated such limits on the Trojan masses in the four 
eccentric systems that we analyzed in this work. We find 
the mass limits to be zero for HAT-P-2b and HD 17156b 
(in agreement with the results of our 3-body integrations that 
were described in §12), 105 M ffi for XO-3b, and 3030 M e 
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TABLE 4 

Observational Constraints on Trojan Masses 



System 


At 




Upper bound (95.4 % confidence) 




[min] 


[M ffi ] 













CoRoT-Exo-1 

CoRoT-Exo-2 

GJ436 

HAT-P-1 

HAT-P-2 

HAT-P-3 

HAT-P-4 

HAT-P-5 

HAT-P-6 

HAT-P-7 

HD 149026 

HD 17156 

HD 189733 

HD 209458 

TrES-1 

TrES-2 

TrES-3 

TrES-4 

WASP-1 

WASP-2 

WASP-3 

WASP-4 

WASP-5 

XO-2 

XO-3 



+0020 
+0012 
-0005 
+0126 
-0055 
-0139 
-0012 
-0025 
-0037 
-0000 
+0026 
-0018 
-0008 
+0002 
-0004 
-0008 
-0014 
-0125 
-0017 
-0122 
+0016 
-0004 
-0014 
+0034 
-0057 



£+0042.2 

oME 

°-0015.7 
c+0034.1 
J -0031.0 
9+0069.6 
z -0082.7 
n+0031.7 
u -0031.7 
9+0138.3 

5 mi 



-0171.0 



"-0009.5 
n+0028.6 

;« 

-0084.6 
e+0010.7 

"+8818:2 

•'-0009.3 
4+0012.7 
-0011.2 
£+0015.2 

dm 

u -0080.9 
£+0054.5 
°-0049.5 
n+0063.8 
'-0052.7 
1+0011.7 
1 -0012.9 
7+0014.5 

\mz 

s+ffl:? 

n+88^8:§ 
-0067.2 



+0018 
+0039 
-0000 
+0022 
-0146 
-0045 
-0004 
-0011 
-0012 
-0000 
+0004 
-0004 
-0007 
+0000 
-0001 
-0006 
-0034 
-0049 
-0007 
-0082 
+0023 
-0006 
-0021 
+0009 
-0384 



n +0061.4 

:mu 

"-0048.0 
4+0001.8 
^-0001.2 
1+0021.0 

-0014.3 
/+0086.9 

-0079.0 
7+0046.4 

4« 

(MH 

°-0097.0 
O+0014.2 
°-00I2.9 
7+0007.5 
-0006.6 
g+0024.6 

7+88§8:l 
6+8883:! 

°-0002.7 
£+0004.6 
"-0004.6 
2+0010.6 



J -0057.5 
q+0031.1 

-0034.4 
O+0042.I 
U -0070.I 
q+0019.3 

-0017.5 
c+0020.0 

Mi 



< 0116.5 

< 0153.0 

< 0002.8 

< 0052.9 

< 0280.5 a 

< 0262.8 

< 0016.5 

< 0147.7 

< 0198.9 

< 0026.5 

< 0016.6 

< 0043.8 a 

< 0022.1 

< 0006.1 

< 0009.9 

< 0024.8 

< 0081.3 

< 0143.8 

< 0070.5 

< 0199.1 

< 0056.1 

< 0043.0 

< 0054.7 

< 0088.2 

< 1097.6" 



<0.35 
<0.14 
<0.12 
<0.32 
<0.10 

< 1.25 

< 0.08 
<0.46 
<0.59 
<0.05 
<0.14 
<0.04 
<0.06 
<0.03 
<0.04 
<0.06 
<0.14 
<0.49 
<0.24 
<0.74 
<0.10 

< 0.11 

< 0.11 
<0.49 
<0.26 



a For this system the upper limit due to dynamical stability is more constraining than that obtained from the data analysis in this paper using the FG method. 
Dynamical stability constraints do not allow for Trojan companions to HAT-P-2b and HD 17156b. And, for XO-3b, the upper limit due to stability is 105 Mjg. 
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FIG. 3. — 95.4%-confidence upper limits on masses of Trojan companions. 
The systems are ordered from the weakest to the strongest upper bound, from 
left to right. 



FIG. 4. — 95.4%-confidence upper limits on Trojan-to-planet mass ratios. 
The systems are ordered from the weakest to the strongest upper bound, from 
left to right. 



for GJ 436b. Thus, for HAT-P-2b, HD 17156b, and XO-3b, 
the upper limit on mj based on considerations of dynamical 
stability is more constraining than the empirical upper limit 
using the FG method. For GJ 436b, the upper limit on mj us- 
ing the FG method ( 2.8 M®) is much stronger than the upper 
limit imposed by the stability requirement. 

5. DISCUSSION 

Exoplanetary science has provided enough surprises that an 
appropriate maxim for observers is: If you can look for a 
novel effect or phenomenon that is at least physically plau- 



sible, then you should do so, especially when this can be done 
with existing data. We have obeyed this maxim by conduct- 
ing a search for Trojan companions to 25 transiting planets. 
Specifically we have put the technique of FG into practice 
with a much larger ensemble than has been previously ana- 
lyzed. We have conducted a search for planets in particular lo- 
cations (L4/L5) with a median sensitivity of ~56 M® , without 
gathering any new data. Instead, we asked the RV data: when 
should the transit occur if there is no Trojan companion? Then 
we consulted the photometric ephemeris to determine when a 
transit actually did occur, and interpreted the time difference 
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as a measurement or constraint on Trojan companions. Our 
results must be understood as constraints on the imbalance 
of mass residing at the L4 and L5 positions. Equally massive 
Trojan companions at those positions would produce opposite 
effects, and no net FG signal, when averaged over the libra- 
tion periods. 

For some systems such as HAT-P-3 and WASP-2, the ex- 
isting RV data are sparse and noisy enough to allow only the 
barest constraints on Trojan companions, with masses com- 
parable to the planetary mass. These constraints are never- 
theless physically meaningful, in the sense that Laughlin & 
Chambers (2002) have shown that equal-mass planets in a 1:1 
mean-motion resonance and circular orbits can be dynami- 
cally stable. In one case, WASP-2, we found a near - 2a ev- 
idence for a timing offset that could be interpreted as a Tro- 
jan companion. This is not compelling evidence, especially 
given that we examined a total of 25 systems, but this system 
is worthy of follow-up. In one case, GJ 436, we have found 
a 2a upper limit of 2.8 M ffi on nip. A positive detection at 
this level would have represented the least-massive planet de- 
tection to date, which is remarkable considering that we did 
not gather any new data. In no case was there evidence for a 
timing offset at the 2a le vel. 

As explained in § 13.31 as part of this study we assessed the 
justification for assuming that a given planetary orbit is cir- 
cular, given the existing RV data and reasonable estimates 
of the stellar age and the timescale for tidal circularization. 
One part of this assessment was the determination of empir- 
ical constraints on ecosw and esinw based on the RV data. 
These results may be interesting to other investigators, inde- 
pendently of our results on Trojan companions. The compila- 
tion in Table 2 of the results for r + and r c ; rc may be useful to 
those who are interested in making inferences about the tidal 
circularization process from the ensemble of transiting plan- 
ets (see, e.g., Rasio et al. 1996, Trilling 2000, Dobbs-Dixon 
et al. 2004, Jackson et al. 2008, Mazeh 2008). The limits 
on ecoso; are also useful for bounding the possible error in 
the predicted times of occultations (secondary eclipses), us- 
ing the relationship given in Eqn. dill . In addition, although 
the planetary radius that is determined from transit photom- 
etry depends mainly on the observed transit depth, there is 
a secondary dependence on the transit timescales (the total 
duration, and the duration of ingress or egress) and the sky- 
projected orbital speed of the planet during the transit. The 
latter quantity is not directly observable; it depends on the 
orbital period, the stellar mass, and the orbital eccentricity 
and argument of pericenter. Thus there is a secondary de- 
pendence of the inferred planetary radius on e and ui (see, 
e.g., Barnes 2007, or McCullough et al. 2008 for a particular 
example). The results of Table 2 can be used to bound the 
systematic error that could arise from this effect. 

Having completed this survey using the existing transit 
data, one may wonder about the achievable limits on Trojan 
companions using data from ambitious future transit surveys. 
We consider here the particular case of the Kepler satellite 
mission (Borucki et al. 2008), whose primary goal is the de- 
tection of Earth-like planets in the habitable zones of Sun-like 
stars. Kepler is also likely to find larger planets such as gas 
giants in the habitable zones of their parent stars, and such 
planets may have lower-mass Trojan companions that are per- 
haps "more habitable" than the gas giants. Such companions 
might be detectable photometrically if they are very nearly 
coplanar with the transiting planet. However, even if they do 
not transit, they can be detected with the FG method. It is 



therefore natural to ask what constraints on Trojan compan- 
ions will be possible for a given planet that Kepler detects in 
the habitable zone of a Sun-like star, using only the photo- 
metric data and the RV data that are routinely gathered for the 
purposes of confirming and characterizing transiting planets. 



m„ = 100 M E=rlh 
m, < 26 M Eo „„ 
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FIG. 5. — Simulated radial velocities and the corresponding 2a upper limits 
on Trojan companions to potential Kepler detections in the habitable zone 
around a Sun-like star with V = 12. 



To answer this question, we consider four different cases, 
in which Kepler finds a planet of mass 300 M® (case A), 
100 M® (B), 30 M® (C) or 10 M® (D), orbiting a star of solar 
mass and apparent magnitude V = 1 2, with a period of 1 yr and 
an orbital eccentricity of 0.1. We simulate RV data for each 
system with a v = 1 m s" 1 , and a number of data points N v that 
seems realistic for the Kepler follow-up program. For cases 
A and B we assume N v = 15. For case C we assume N v = 30, 
which is sufficient to measure the planetary mass to at least 
10% according to the expression for the signal-to-noise ratio 
from Gaudi & Winn (2007), 



S/JV-0.6 



nip 

Mo, 



100 



1/2 



-1/2 



m s 



-1/2 



10 



-0.2(V-12) 



D 

3.6m' 



(16) 



which was intended to approximate the case of the HARPS 
instrument on the ESO 3.6 m telescope. For the challenging 
case D we assume N v = 70, corresponding to a 20% uncer- 
tainty in the planetary mass. We perform a MCMC analysis 
on the simulated data for each system, just as was done for the 
25 transiting systems in this study, and obtain the correspond- 
ing constraints on Trojan companion masses. 

Fig-Hshows the simulated data and the results for the four 
cases. We find the 2er upper limits on the mass of Trojan 
companions to be 25 M®, 26 M®, 21 M®, and 24 M®, for 
cases A, B, C and D respectively. The results are limited by 
the RV follow-up program; the superb photometric precision 



Trojan Companions to Transiting Exoplanets 



9 



of Kepler does not lead to correspondingly superb constraints 
on Trojan masses. 

The results for cases A-D are all of the same order of mag- 
nitude. This is because Trojan detectability depends primar- 
ily on <x v / \/N\,, which only varies by a factor of 2.2 between 
case A and D. For small mp/mp and large £ (see Eq. 9), 
Trojan detectability is indeed independent of planet mass for 
fixed <7 V / \/N~ v . A larger planet produces a larger RV semi- 
amplitude, and hence offers greater sensitivity in measuring 
Af, but the conversion from Af to mp varies inversely with 
planetary mass. 

However, for cases C and D, this scaling is not precisely 
obeyed. One reason is that for the larger values of mp / mp 
that are relevant in those cases, the relation between mp and 
Af is nonlinear (see Eq. [2]), and therefore the error in mj is 
not Gaussian. For case D, there is an additional source of 
non-Gaussianity: the signal-to-noise ratio is low enough that 
the correlations between Af, K, e, and u> become important. 
This means that the upper limit on nip is less constraining than 
one would predict based only on <r v / \/r%. 

We conclude that a "serendipitous" search for Trojan com- 
panions to the habitable-zone planets that will be detected and 
confirmed as part of the Kepler program will be sensitive to 
planets of approximately Neptunian mass or larger. Of course 
the sensitivity could be improved by obtaining additional RV 
data as part of a more focused search effort. 

Considering the HARPS spectrograph (Pepe et al. 2002), 
mounted on the 3.6 m ESO telescope, as a fiducial instru- 
ment for precise RV measurements, and assuming the noise 
to be limited by photon-counting statistics, the measurement 
uncertainty can be obtained by scaling current results (Lovis 
et al. 2005, Gaudi & Winn, 2007): 



where V is the apparent visual magnitude of the star, and D is 
the aperture of the telescope. Here we have assumed a 60 min 
exposure and a G-type star. 

For small Trojan-to-planet mass ratios, and assuming a 
circular orbit, one can determine a nominal estimate of the 
sensitivity of current observational facilities to detect Trojan 
companions to Kepler planets. Under these asumptions, the 
signal-to-noise ratio At /a a, for detecting a Trojan compan- 
ion is given by: 



S/JV~0.56 



m s 



M e A 100/ \AUJ 

1/2 

3.6m 



-1/2 
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-0.207-12) 



(18) 



1q0.2(V-12) 

= „ , — m/s, 



D/3.6m 



(17) 



Let us assume a Kepler detection of a Jupiter-mass transiting 
planet at 1 AU around a sun-like star with V = 12, observed 
with the HARPS instrument. Then, considering 100 obser- 
vations evenly spaced in orbital phase, a 3er Trojan detection 
(i.e. Af > 3 a At) can be made, if a Trojan mass imbalance 
°f m T 5.36M e existed in the orbit. If the planet were or- 
biting instead at 0.03 AU (a "hot Jupiter" orbit), then a de- 
tection would be possible for mp > O.93M0. It would seem 
that searching for Trojan companions is a promising alternate 
channel for finding small and potentially habitable bodies in 
the habitable zones of their parent stars, even if the transiting 
planets themselves are too massive to be habitable. 

We thank Jack Wisdom, Scott Gaudi and Eric Ford for help- 
ful conversations. We further thank Scott Gaudi for providing 
a detailed and helpful review of the manuscript. We are grate- 
ful to the William S. Edgerly Innovation Fund for partial sup- 
port of this work. 
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